The incommensurate magnetic structure of YMn205: a stringent test of the 

multiferroic mechanism 
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We have determined the magnetic structure of the low-temperature incommensurate phase of mul- 
tiferroic YMn205 using single-crystal neutron diffraction. By employing corepresentation analysis, 
we have ensured full compliance with both symmetry and physical constraints, so that the electri- 
cal polarization must lie along the b axis, as observed. The evolution of the spin components and 
propagation through the commensurate-incommensurate phase boundary points unambiguously at 
the exchange-striction mechanism as the primary driving force for ferroelectricity. 
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The family of compounds with general formula 
RMn2 05 {R=Y, Rare earth, Bi and La) [ij occupies a 
special place among the so-called "novel" multiferroic 
materials, in which ferroelectricity is directly induced 
by the onset of inversion-breaking magnetic ordering. 
In fact, whereas for most other multiferroics, such as 
TbMnOa 0, Q Ni3V2 08 4] and many other materi- 
als, the so-called spin-orbit (SO) mechanism 0, Q has 
clearly been identified as the primary driving force for 
the electrical polarization, the origin of ferroelectricity 
in RMn2 05 remains controversial. Three aspects of the 
RMn2 05 phenomenology suggest that the SO mechanism 
may not play a major role : the direction of the electri- 
cal polarization cannot be rotated by an applied mag- 
netic field 0]; the main ferroelectric phase is magnet- 
ically commensurate while the low-temperature incom- 
mensurate phase (LT-ICP) is only weakly ferroelectric, 
and, most importantly, in the ferroelectric commensurate 
phase (CP), moments in the a&-plane are almost collinear 
(non-coUinearity is a strict requirement of the SO mech- 
anism). An alternative explanation, based on the simple 
exchange-striction (ES) mechanism in the context of an 
acentric quasi-coUinear structure, has been shown to ac- 
count qualitatively for the ferroelectric behavior of the 
different magnetic phases [1, and recent electronic 
structure calculations not including the SO interaction 
[iQI] have reproduced the observed value of the electri- 
cal polarization P for certain values of the parameters. 
However, up to now the SO mechanism could not be com- 
pletely ruled out, since accurate neutron diffraction de- 
terminations of the CP magnetic structure [ll|, [l2| have 
identified a small cycloidal component in the ferroelectric 
phases. 

One striking experimental fact about RMn2 05 is the 



sudden, dramatic reduction in the value of P at the CP 
to LT-ICP transition. Understanding how the different 
components of the magnetic structure change could pro- 
vide key insight into the multiferroic mechanism. More 
specifically, we have shown [l3| that, starting from the ex- 
perimental magnetic structure, one can construct time- 
reversal-even polar vectors .specific to each mechanism, 
which are directly proportional to the SO and ES contri- 
butions to P. Most likely, the coupling constants do not 
change across the transition; therefore, the magnitude of 
the electrically active polar vector must, like P, be greatly 
reduced across the transition. For this approach to be 
valid, the magnetic structures must be known with great 
accuracy. In this respect, YMn2 05 is an ideal model 
system since there is no rare-earth magnetic ordering in- 
volved. At present, whereas the YMn205 CP magnetic 
structure is well established, and has been validated inde- 
pendently by two groups [111, [l^ 7 the LT-ICP structure 
is not known with confidence. The two structures pre- 
sented to date from neutron powder diffraction [l3 |. and 
more recently, from single crystal data [15[ are signifi- 
cantly different, and, crucially, neither of them exploits 
the full set of symmetry and physical constraints, requir- 
ing that the residual electrical polarization in the LT-ICP 
remains parallel to the b axis. 

Here, we present a single crystal neutron diffraction de- 
termination of the YMn205 incommensurate magnetic 
structure that is fully compliant with symmetry and 
physical constraints. In particular, we have employed 
corepresentation analysis to impose the ■2- point-group 
symmetry, which is lower than the CP symmetry (m2m) 
but still guarantees that all the contributions to P are 
directed along the b axis, as observed experimentally. 
The fit to the data, collected on two independent mag- 
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netic domains, is much better that for previous structures 
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15j , and is essentially equivalent to our best uncon- 
strained refinement. The spins responsible for the ES 
polarization are now found to be nearly perpendicular 
and are modulated with opposite phases, leading to a 
reduction of the ES polar vector by a factor of 9 with 
respect to the CP. On the contrary, the magnitudes of 
the SO polar vectors increase significantly in the LT-ICP, 
while a new polar vector of SO origin and of smaller mag- 
nitude develops due to the formation of long-period cy- 
cloids along the a axis. Barring an accidental cancelation 
of SO terms, which is extremely unlikely to occur for all 
RMn205 compounds and all temperatures, this proves 
conclusively that ferroelectricity in the RMn2 05 com- 
mensurate phases is predominantly of exchange-strictive 
origin. However, it is possible that the SO mechanism 
may provide a contribution to the residual polarization 
in the low-temperature incommensurate phase. 

The YMn2 05 single crystal for this study is the same 
as for reference , and the growth protocol is described 
therein. The experimental setup was also the same but 
the data were collected at 2K. The main difference be- 
tween the experiments on the CP and LT-ICP phases is 
at the data reduction stage: in the CP with propagation 
vector k — 27r(i,0, |;), there are two inversion-related 
domain that scatter onto the same points in recipro- 
cal space, yielding single Bragg peaks on the position- 
sensitive detector (PSD). On the other hand, the LT- 
ICP has 4 distinct domains related by inversion and/or 
rotation around the c axis. The rotation-related do- 
mains have propagation vectors fci = 27r(0.48, 0, 0.28) 
and k2 = 27r(-0.48, 0, 0.28). Scattering from these do- 
mains occurs at close positions in reciprocal space, yield- 
ing split peaks on the PSD. An ad hoc piece of software 
was developed to de-convolute the peaks from the two do- 
mains and exclude the reflections that were completely 
overlapped. Because of the different integration meth- 
ods employed for the CP and LT-ICP, absolute scaling 
of the magnetic moments was verified using the neutron 
powder diffraction data reported in [l4'|. The FuUProf 
program [l6l | was used for the refinements of the mag- 
netic structure. In the final analysis, reflections from the 
two domains were combined in a single refinement, with 
the spin components of the two domains constrained by 
symmetry (see below). 

For the corep analysis, we use the Kovalev conventions 
as in lljj. The propagation vector is /c = 27r(0.48, 0, 0.28), 
(/c3 in the Kovalev notation). The little irrep group con- 
tains two operators — the identity hi and the a-glide /i27 
perpendicular to h. According to Kovalev (see Table [l|, 
there are only two irreps, both 1-dimensional, each gener- 
ating a corep. Both sites, Mn'^+ and Mn'*+, split into two 
orbits in the little group, which are then recombined by 
the coreps. The transformation matrix P is the identity. 
The axial-vector modes were symmetrized with respect 
to the center of symmetry located at position 0, 0.5, 0.5. 



The two-fold screw h^ is located at 0.5,?/, 0.5. In this 
respect we depart from Kovalev's conventions. Magnetic 



TABLE I: Small irreps (A) and coreps (D) of space group 
Pbam for propagation vector — (/i, 0,/i). The symme- 
try operators are in the Kovalev notation, e — e~^'^* ^ — 
0.0628 - i 0.998). 

hi /i27 Kh25 Khs 



Ai/Di 
A2/D2 



structures built from single- corep modes are by construc- 
tion invariant by inversion (/125), and are also invariant 
by rotation around the b axis (/13) and reflection through 
a plane perpendicular to the 6-axis (/127), except for the 
phase factor e or — e, which, in an incommensurate struc- 
ture, is always equivalent to a translation and does not 
affect the point-group symmetry. Our aim is to construct 
a magnetic structure that allows the development of po- 
larization along the b axis — we want to drop the /125 
and /i27 invariance but retain the h^ rotation around b 
(point group symmetry •2-). Inspection of Table [T] imme- 
diately suggest that this is accomplished by the combi- 
nation Di + iD2. In fact, inversion and rotation appear 
in the little corep group as antiunitary operators. The 
imaginary unit in front of D2 yields a sign change for 
these operators, ensuring the correct symmetry behav- 
ior. 

As usual, each corep mode for a single spin component 
has two parameters — an amplitude and a phase, the 
latter describing the phase difference between inversion- 
related sites. Therefore, the most generic magnetic struc- 
ture invariant by •2- symmetry is described by 24 param- 
eters (2 parameters x 3 spin components x 2 coreps x 
2 sites). However, this number can be further reduced 
by imposing physical constraints. Sites related by in- 
version in the paramagnetic phase have almost identi- 
cal magnetic environments (except for the tiny ferroelec- 
tric displacements). Hence, we expect not only that the 
wave-amplitudes will be identical, but also that the spins 
follow the same propagation minus a phase factor. This 
reduces the number of phases per site from 6 to 1 and the 
total number of parameters to 14. One of these param- 
eters represents the overall phase and can also be fixed 
(we have chosen to fix ub^. — 0, see below), yielding 13 
parameters in total, a drastic reduction from 48 for an 
unconstrained refinement [l5| . 

The j component of the magnetic moment on site I and 
unit cell n (at a position i?„ from the origin) is written 
as 

m\{n) = S]e'^-^- +C.C. (1) 
The Fourier coefficients have been parameterized as, for 
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example for component x of site h: 

Sl = \{uhl + ^ubl)e^^ 



(2) 



where $ is a global phase common to all components. 
The refined parameters for Domain 1 and the appropri- 
ate transformation to obtain those of Domain 2 are listed 
in Table |TT1 The fit is far superior to that obtained using 



TABLE II: YMn205 LT-ICP magnetic structure at T = 2 
K. Parameterized Fourier coefficients for Domain 1, fci = 
27r(0.48, 0,0.28). The refined values are: uhl = -2.45(3), 
uhi = 0.30(6), uhl = -0.27(10), uhl = 0, uhl = 1-72(1), 
uhl = 1.03(3), ipb = -0.36(1), uai = 3.50(3), uaj = 
0.46(8), ua\ = 0.26(7), ual = -0.53(8), ual = -2.00(1), 
ual = 1.20(4), = -1.650(7). Rf = 6.75/6.67 % and 
7?!' = 10.7/9.83 % for Domain 1/Domain 2, respectively. 
The transformations to obtain the Domain 2 parameters are: 



k2 

,2 



27r(-0.48, 0,0.28) ul 



u^ —* — m;, Uy — > — ui — * ui for both a and h sites, and 
i^b ^ tpb, i^a —> —tpa- Labeling of atoms is as in [Tl| . Sup- 
plementary Table I in the EPAPS contains a complete list of 
parameters for both domains, using the FuUProf conventions. 
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the model by Kim et al. [l5| {Rp =i 6.7 % vs 15.5 % for 
Kim), and is better than any unconstrained refinement 
we could obtain based on a simulated annealing struc- 
ture. The magnetic structure of the LT-ICM is depicted 
in Fig[TJ The top panel shows the a6-plane projection of 
the envelope of the magnetic structure, i.e., several spin 
orientations describing a complete period of propagation 
are shown on the same site. The bottom panel shows 
a perspective view of the magnetic structure illustrating 
the phase relations between spins on different sites. The 
LT-ICM adopts the most generic type of single-A: struc- 
ture, in which spins rotate in an arbitrary plane with an 
elliptical envelope — in other words all spin components 
are involved in the cycloids. These envelopes are strictly 
related by symmetry on the different Mn'^+ and Mn^+ 
sites, and differ only by a phase factor. Mn^+ and Mn^+ 
envelopes are almost co-planar, but the Mn'^+ ellipse is 
much more eccentric (the semiaxes are 3.65 and 2.21 /is 
for Mn3+ and 2.48 and 2.00 /is for Mn''+ ) , refiecting 
the different anisotropy of the two sites. The Mn^+ long 
semiaxis coincides with the pyramid axis, consistent with 
the easy magnetic axis from susceptibility measurements 
and also with the results on the CP, whereas the short 



semiaxis points towards one of the oxygen atoms in the 
pyramid. The most important aspect of this structure is 
the out of phase relation between zig-zag chains running 
along the a axis (bottom panel of Fig. [Ij . This has two 
consequences: first, neighboring spins on different chains 
are essentially orthogonal; second, even when the angle 
is not 90° due to the eccentricity, the dot product across 
chains cancels out almost exactly due to the sign change 
in different part of the structure. As the ES polarization 
is proportional to this dot product, we argue that this 
loss of phase coherence is the primary cause of the loss 
of ferroelectricity in the LT-ICP (see below) . As in many 
other cycloidal magnets, the CP to LT-ICP transition 
can be well explained in terms of competition between 
anisotropy (which favors quasi-collinear arrangements) 
and entropy (which favors large-moment cycloids at low 
temperatures), but, intriguingly, in RMn205 the ferro- 
electric properties of the two phases are reversed with 
respect to TbMnOs [5]. In order to calculate the dif- 




FIG. 1: (color online) Low-temperature incommensurate 
magnetic structure of YMn205 . Top: projection of the mag- 
netic structure onto the afo-plane. Mn'^"'" and Mn*"'" are shown 
with red and green color, respectively. The ellipsoidal con- 
tours, more elongated for Mn''"'" , show the envelop of the mag- 
netic modulations. Adjacent AFM zig-zag chains, labelled I 
and II, are shown by thick black lines. Bottom: Projection 
of the magnetic structure onto the fcc-plane. The black arrows 
represent the normal direction to the ellipoidal modulations 
in chain I and II. 

ferent contributions to the polarization, it is useful to 
extract modules (m^) and phases (0") for the complex 
arrays Za and z^, defined as Zax = ua^ + i ual ^i'^-' 



arg(za) ; 06 = arg{zb) 
\za\ ; rub = \zb\ 



(3) 



Using these equations, one can readily calculate all the 
macroscopic quantities of interest, and in particular the 
polar vectors corresponding to the different mechanisms 
|13| . Since the point-group symmetry is ■2-, only the y 
component of any polar vector will be non-zero. This can 
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readily be verified by calculating the vector component 
explicitly, but it is not shown here. By employing the 
transformations in Tab. |TT1 we can also verify explicitly 
that all the polar vector change sign in Domain 2. The 
formulas given here below are valid for the LT-ICM (no 
umklapp term while for the CP we use eq. 3 and 4 
in The ES polar vector is calculated as: 



= ^'^ajmhj [AcosipbSm^pa sin(0;^^- - (j)'^^)] 



(4) 



The two spin-orbit components that were already present 
in the CP, originating from 6c-plane cycloids and relating 
to SO interactions across the Mn^+ planes {Sy) and Y 
planes (Sy) are: 



sl - 

Sy = mbj^mb^ [2 sin(2'0b + gz) sin( 



%-<f>L)] (5) 



In order to evaluate the spin-orbit contribution due to 
the in-plane projection of the cycloids (which was absent 
in the CP), we only consider the Mn'^+-Mn'*+ within the 
chains, since the other bonds have a small projection on 
the direction of propagation (spins on Mn^^ atoms re- 
lated by inversion are also almost antiparallel). 

Ly = — 4sin(g2:/2 — ipa) cos^pi, cos a • 
[mbxrriay sin((/)^'^ + Ky) + mbyma^ sin(0^^ + (/)^^)] (6) 

Table IIIII lists the polar-vector components for the CP 
ll( and LT-ICP, as determined from the previous for- 
mulas. The relative signs of the different components 
are uniquely established, but the overall sign depends on 
the choice of domain. Tab. IIIII provides an essentially 



TABLE III: Exchange-striction (Ey) and spin-orbit (S^, Sy 
and Ly) polar vectors for the CP and LT-ICP phases of 
YMn205 . All the values are in fi%. Statistical errors are 
propagated from the magnetic structure refinements. 



Phase 


Ey 


sl 


Sy 


Ly 


CP 
LT-ICP 


22.6(5) 
-2.5(3) 


-0.16(4) 
-1.01(8) 


-0.40(5) 
-1.3(1) 




0.7(3) 



unambiguous test for the multiferroic mechanisms: the 
drop of the ES component Ey is by far the strongest can- 
didate to account for the sudden loss of ferroelectricity 
at the CP/LT-ICP transition, as the SO components Sy 
and Sy significantly increase in magnitude. The only 
plausible SO-based alternative is an accidental cancela- 
tion between Sy, Sy and Ly. However, this would be an 
extraordinary coincidence, since the corresponding SO 
coupling constants may not even be of the same order of 
magnitude, as the exchange pathways are completely dif- 
ferent. Furthermore, this coincidence would be required 
for all the RMn205 compounds and all temperatures, in 



spite of the significant perturbations introduced by the 
rare earth magnetism [18| . The origin of the residual fer- 
roelectricity in the LT-ICP is, however, a different mat- 
ter. The drop in Ey (~ 89%) is somewhat larger that 
the drop in P (~ 80%), making it plausible that the SO 
mechanisms may contribute significantly to the LT-ICP 
electrical polarization. 

In summary, we have determined the magnetic struc- 
ture of the low-temperature incommensurate phase of 
YMn205 from single-crystal neutron diffraction data. 
By employing corepresentation analysis, we have im- 
posed strict symmetry and physical constraints, so that 
all component of the electrical polarization lie along 
the 6-axis, as observed experimentally. Quantitative 
analysis of the magnetic structure change across the 
commensurate-incommensurate transition points unam- 
biguously to exchange-striction as the primary origin of 
ferroelectricity in this system. 
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